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Abstract 

We compare ground state properties of 3D Ising Spin Glasses with Gaussian 
couplings with results from off-equilibrium numerical simulations at non zero 
(but low) temperatures. We find that the non-zero temperature properties of 
the system smoothly connect to the T = behavior, confirming the point of 
view that results established at T = typically also give relevant information 
about the physics of the system. 



1 Introduction 



Spin glass physics is difficult. Because of how complex it is to establish results that go 
further than the mean field level (see for example |l| and references therein) frequently 
numerical simulations are the technique of choice (see for example || and references 
therein) . 

Recently ground state techniques have become popular (see for example [[| and 
references therein). The big advantage here is that one is (down in the broken phase) 
surely as far as possible from the critical temperature T c . Our goal is to observe the 
effects of the T = fixed point, and to remove spurious effects due to the critical 
point at T = T c : these contaminations are minimized when working at T = 0. Also 
the use of exact ground states solves the problem of thermalization, that is typically 
very severe on spin glasses at low values of T. The price one pays is that contacts 
with T ^ physics are not clear: it could even be (but it does not seem to be so, 
and here we will help in showing that this is not the case) that what happens at 
T = is essentially different from what happens for an even infinitesimal value of the 
temperature. 

Here we will use a further tool to check consistency of the T = and the T ^ 
phase space. Thanks to off-equilibrium dynamical simulations we will be able to 
study large lattice sizes and restricted sectors of the phase space: this will allow to 
show better that T = physics is compatible with the results one obtains at finite 
values of the temperature. 

The second important part of this work is that, always thanks to off-equilibrium 
dynamical numerical simulations, we are able to keep under control crucial observ- 
ables like overlap-overlap correlation functions and block overlaps. These observables 
are very important since they allow to discriminate potentially misleading situations 
like the ones where interfaces induce a seemingly non-trivial overlap probability dis- 
tribution P(q) from a true replica symmetry breaking. 

2 Model, Algorithm and Observables 

We have simulated a three dimensional (3-D) Ising Spin Glass with Gaussian couplings 
on a cubic lattice with periodic boundary conditions. The Hamiltonian is 



where the sum runs over all couples of nearest neighbors, the o".j = ±1 are Ising spins 
and the couplings Jjj are quenched random Gaussian variables with zero mean and 
unit variance. 

We have mainly focused on the measurements of two very important observables: 
they are both useful to distinguish situations where a true non-trivial behavior of 
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the probability distribution of the overlaps is present from situations where interfaces 
could present a misleading situation (because of finite size effects or even in the infinite 
volume limit, see the discussion in |4], [j| and references therein) 

The first observable is an overlap, i.e. the measurements of how similar two typical 
configurations at equilibrium are. Here we define the overlap only in a small part of 
the lattice: we compute the block overlap qg on a 2 x 2 x 2 cube (that we call a 
£>-cube, here with B — 2): 




ieB 

with the usual site overlap 

Qi = (3) 

where a and r are two independent equilibrium configurations of the system in the 
same realization of the quenched disorder. We call, as usual, q the total overlap 
computed on all the lattice, i.e. the average of the site overlap qi over all lattice sites. 
We notice also that in the following we will, as usual, denote with brackets the thermal 
average for one given realization of the quenched disorder, (( ■ )): when averaging at 
a given Monte Carlo time t we will mean averaging over different realizations of the 
dynamical process, and we will indicate this average with (( • )) t . We will denote 
with an over-line the average over the quenched disorder, ( ■ ). 

We have computed the full probability distribution Pg(gH, t, T) of the B = 2 block 
overlap q& as a function of the temperature T and the (Monte Carlo) time t of a 
so-called off- equilibrium dynamics. In an off-equilibrium numerical simulation (see 
for example ||) we use a "very large" lattice size, and we work at values of the 
temperature low enough to make sure that the system is not thermalized: this is 
true already for medium size lattices in the case of system that are characterized by 
slow dynamics, as is the case here. In this situation we work by extrapolating off- 
equilibrium, finite time measurements to infinite time. The advantage of the method 
is on one side in the fact that we are basically free from finite size effects, and on the 
other side that we can learn additional information from the pattern of the approach 
to equilibrium (for example in || this turns out to be a nice way to learn about the 
minimal value of the allowed overlap, q m i n )- On the other side the disadvantage is 
that we have to rely on an extrapolation to infinite time, that can turn out to be not 
so trivial. 

The second relevant observable is the overlap-overlap correlation function (defined 
on the full lattice) computed at distance r, and it is defined as 

C(r,t,T) = { qi q i+r ) t . (4) 

Here we will only analyze C(l,t, T), since among all the different C for different r 
values it is the less affected by statistical errors. 
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In our working conditions C(r,t,T) depends on the initial total overlap prob- 
ability distribution P(q,t = 0,T). In our case, since we have used very large lat- 
tice sizes, since we are starting from an average zero overlap (we start by evolving 
two independent, uncorrelated random configurations) and since we are limited to a 
number of Monte Carlo sweeps (of the full lattice) of the order of 10 7 , we are not 
able to change away from zero the total overlap, and at least roughly we have that 
P(q,t,T) ~ P(q,t = 0, T). Because of that we have that the correlation functions 
we measure are very different from the equilibrium correlation functions, and can be 
written as 



where by C(r,t,T,q) we denote the spatial correlation functions at time t and tem- 
perature T computed by including only couples of configurations having total overlap 
q, and by C (r,t,T,q) we denote the correlation functions computed with initial 
conditions P(q,t = 0, T). 

In the numerical simulations that we will discuss in the following we have been 
starting from an initial probability distribution peaked around zero total overlap, i.e. 
our results are determined by the value of C(r,t,T,q = 0). In other terms we are 
discussing here the dynamics of the q = sector of the system (see || for further 
details). As we have stressed in the introduction, our main goal here is to relate 
numerical results computed from an off-equilibrium dynamics at non zero temperature 
with equilibrium results obtained at zero temperature using techniques that give exact 
ground state configurations on reasonably lattice sizes (see for example 0,0). 

In the recent studies computing spin glass ground states the technique of choice 
was based on computing the so-called mixed overlap, i.e. the overlap computed be- 
tween a ground state in a finite volume with periodic boundary conditions (pbc) in all 
directions and another ground state obtained when imposing anti-periodic boundary 
conditions (abc) in one direction and pbc in the other directions of the cubic lattice 
with the same realization of the quenched couplings J. In a theory with continuously 
broken replica symmetry |9[] in the infinite volume limit there are many ground states 
with total site overlap q and total link overlap q\ (that coincides with the correlation 
function C computed at distance 1) different from one ||. On the contrary a droplet 
approach (TTJ] implies a behavior similar to the one of an usual ferromagnet, where 
the link overlap (and the site overlap if we activate an infinitesimal magnetic field) 
go to one in the thermodynamical limit. Both and || compute the B = 2 block 
overlap. 

Now a key remark is in order. The procedure we have just described to compute 
overlaps in the ground state sector leads to the comparison, if the RSB picture holds, 
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of ground states with q = (changing the boundaries we favor ground states with 
small overlap): that means that these results based on ground states describe the 
link overlap in the q = sector of the system. This is the main reason of interest of 
these dynamical simulations: we are simulating here the q = sector of the theory, 
and since we are working on very large lattice sizes we are only sampling this sector. 
In other words when extrapolating our measurements of 0^(1, T) (by the subscript 
oo we mean we have already extrapolated, hopefully faithfully, to t — oo) from our 
low values of the temperature T down to T = we should find, if the RSB picture is 
consistent, the same results found in the direct analysis of the ground states. 

We will also try to check the behavior of a different relevant observable. In || 
the authors have computed the probability that, for example, a B = 2 block does 
not intersect an interface: one can relate this probability to the probability of finding 
that a B = 2 block is ±1. Again our strategy will be to compute the probability of 
this event at non zero temperatures and to extrapolate the result down to T = 0. 

A very direct relation about the second moment of the (B = 2) block overlap and 
the q — q correlation function |] can be established by noticing that 

{<&) (t, T) = l \C(0, t, T) + 3 C(l, t, T) + 3 C(V2, t, T) + C(V3, t, T)} . (6) 



3 Numerical Results 

Our numerical simulations have been run on a APE100 parallel supercomputer [|Pl| : 
on the more powerful version of the computer (the so-called "tower" version) our 
Monte Carlo code for simulating 3D Ising Spin Glasses with Gaussian couplings has 
a peak performance of 5 GF. 

For all the numerical simulations described in this note we have averaged our 
data on four different realizations of the random quenched disorder. In each disorder 
realization we have simulated two independent copies of the system (useful to compute 
the different overlaps, block overlaps and correlation functions). We have used a 
straightforward Metropolis algorithm. We have used a lattice of linear size L = 64 
and volume V = L 3 = 2 18 . 

We have analyzed two kind of runs. "Short" runs have been used to compute the 
extrapolated value of the overlap-overlap correlation function. In these short runs 
we have simulated temperatures T = 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.35 and 0.3 during 
819 200 Monte Carlo steps. We have presented a preliminary and partial analysis of 
some of these data in reference (only for T > 0.35 and only as far the dependence 
of the critical exponents over temperature is concerned: we had not analyzed, for 
example, the prefactors, which are very relevant to the physics studied in this paper, 
see below). 

In the "short" runs we use a simple scheme. We start from two initial independent 
random configurations, and suddenly start iterating at the working temperature well 
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below the critical one: T < T c (T = 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.35 or 0.3). At each T 
value we will eventually extrapolate to infinite time, to get a T dependent asymptotic 
value, that we will, in turn, try and extrapolate down to T = 0. 

In order to compute the dynamical behavior of the block overlap we have used 
"long" runs. Here we have used the previous scheme for four values of the tempera- 
ture, T = 0.7, 0.6, 0.5 and 0.35, with 6 553 600 Monte Carlo updates of the full lattice 
performed at each value of T. 

In the "short" runs we measure the time dependent observables on logarithmic 
time scales: i.e. at times t = 100, 200, 400, 800, ... In the "long" runs we measure 
every 32768 steps. 

In the next two subsections we will discuss separately our numerical results for 
the correlation functions and for the window overlap. 



3.1 Correlation Functions 

In this section we will always indicate by C(l,t, T) the q = component of the 
overlap-overlap correlation function, i.e. C(r = l,t,T,q = 0), that, as we have 
discussed in detail in section ^|, is the quantity we are measuring. 

The dependence of C(l,t,T) on time can be expressed through the knowledge of 
the dynamical critical exponent z(T) (that depends on T). We have that 



C(l,t,T) = f(t^m,T 



From numerical simulations [[12], |13|, |Tj, [15], |1| (that here are in very nice agreement 
with real experiments []T7|) we know that 



z(T) = - (6.2 ± 0.3) . 
Another useful piece of information Jl4], [H| [16| is that 



(7) 



C(x 1 t,T) oc -^exp 



x 



x 



t{t,T) 



3/2' 



(8) 



where £ (T) oc t 1 ' z ^ T ' is the time dependent dynamical correlation length of the system. 
Because of that we expect that 



logC(l, t, T) ~ a(T) + b(T)t 



-3/(2*(T)) 



(9) 



asymptotically for large time. As t diverges C(l,t,T) tends to the constant value 
exp (a(T)), that we have called prefactor in the previous section. 
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In figure [l] we show the extrapolation of C(l, t) to t = oo for the lowest temperature 
we use, T = 0.3 (this is the most difficult extrapolation, and fits at larger T values 
are easier). The dashed line is for the best fit to the form 

log C(l, t, T) = a(T) + b(T) t~ c ^ . (10) 

The fit is good, and the extrapolation looks reasonably safe. By this kind of fits 
we are able to obtain reliable extrapolated values for all the values of T we study, 
that we show in figure ^| versus the temperature T. We have checked that the value 
of c(T) of the best fit to the form ITU] is in very good agreement (i.e. inside the 



margins suggested by the estimated statistical error) with the value of 3/(2z(T)). 
The discrepancies between the extrapolated value obtained by using a three or two 
parameter fit (i.e. by fixing c(T) to 3/ (2z(T)) in equation pjor by letting it free) are 
inside the statistical errors. 



For a detailed numerical study of the functional form suggested in equation [TTJ we 
refer the reader to references 0, [15], An exhaustive analysis of the correlation 
functions computed at zero total overlap has been performed in three [13, [IH |TB 



four and six dimensions []19| : in all cases the system turns out to have a very 



similar behavior, independently of the spatial dimensionality. In six dimensions the 
agreement with the quantitative analytical predictions obtained by De Dominicis, 
Kondor and Temesvari [^Uj is very good |05] . 



In the same figure we have also plotted the statistical interval allowed for the value 
of the link overlap at T = from ground state calculations (see ||). The statistical 
error on our values of 0^(1, T) grows when T decreases, but in the limit given by 
these errors the T = result and a reasonable extrapolation of the finite T results 
obtained here look completely consistent. 

The interval we have drawn in the figure is 0.79 ± 0.07, and we have computed 
it by looking at data from ||. We have used the four different (but statistically 
compatible) values obtained there for q\. By extrapolating of q\ for L — > oo one 
obtains qi = 0.755 ± 0.015 or qi = 0.80 ± 0.06 using a linear or a quadratic fit 
in 1/L, respectively. By studying correlation function reference || also obtained 
C(l) = 0.732 ± 0.008 (for the transverse correlation) and C(l) = 0.722 ± 0.005 (for 
the perpendicular correlation). Taking into account all these figures we have obtained 
the confidence interval (at one standard deviation) for the link overlap that we were 
quoting before. 



3.2 Block Overlaps 

Our second set of measurements concerns the so called block overlap (that we also 
call window overlap or box overlap, and have been computed onaSxfixB = 
2x2x2 box). As in the case of correlation functions we start from two independent 
random configurations, and cool them down according to the annealing schedule we 
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Figure 1: Values of C(l,t,T) at T = 0.3 (the lowest value of T that we use) versus 
£-3/(22) _ a j so s } low our b es t fit t the form ITOl 



have described before (using what we have called the "long" runs). During these runs 
we have computed the full probability distribution of the window overlap at different 
times t. 

In this section of this note our main goal is to compute Pg(±l), i.e. the value 
(extrapolated to infinite time first and than T = later) of the probability that a B 
block has overlap of maximal modulus. To simplify the notation from now on we will 
be working with a symmetrized overlap probability distribution, i.e. we will consider 
Pg(|g|,t, T), that can have support in (0,1): this object is interesting because it is 
exactly the probability that the interface between different phases does not intersect 
with a cube of size 2. As we have already discussed this probability has been computed 
in reference |§ at T = 0, and we will try here to connect these results with the T^0 
physics. 

The main result of reference |§ is that, at T = and in the limit of the lattice 
size L — > oo, 

Pi {K) ~ 0.65 ±0.05 . (11) 

Our working Ansatz for the time dependence of Pb(|?| = 1) has been, as we will 
justify now, that 

logP B (|g| = l,t,T) =logPr°°(|g| = 1, T) + -j^JL , (12) 
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Figure 2: Values extrapolated to infinite time of the correlation function at distance 
d = 1, Coo(l,T), versus T. We have also marked by two horizontal dotted lines 
the interval where the value computed in reference || using T = ground state 
calculations lies. The consistency of the two results is clear. 



and we have used z{T) = ^r. 

Let us see why things work in this way. From equation || we know that < q 2 >b 
scales as an overlap-overlap correlation function. But for large times, and for distances 
such that x -C t x l z equation |8| tells us that the correlation function behaves as 



C(x, t) oc —= 



A 



x 



t 3/(2z(T)) 



where A is a suitable constant. This implies in turn that 



<q 2 B >(t,T)~a(T) + 



b(T) 



t 3/(2z(T)) ' 

where a(T) and b(T) are constants that only depend on T. Now since 



<q 2 B >(t,T) = j ^dqP B (q,t,T) q 2 



if we define AP B (q, t, T) = P B (q, t, T) - P B =°°(q, T) it is clear that 



dqAP B (q,t,T) q 2 



t 3/(2z(T)) 



(13) 



(14) 



(15) 



(16) 



Things are more clear if we write equation [16] for discrete values of q& (that is what 
happens for finite £>): 

J2 AP B fa, t, T) q 2 + AP B (q = 1, t, T) = . (17) 

<7<1 

If we assume that all the components with q > scale in the same way this shows 
that 

AP B (q = 1, t, T) = P B (q = 1, t, T) - P B =°°(q = 1, T) ~ r 3 ^ 2 *™ . (18) 

and so, 

P B (|g| = l,t,T) = P^Qql = 1,T) + , (19) 



and we recover [12| taking logs in [19|. We have used [12| instead of [19] since we have 
(slightly) better scaling even for small times. We have checked that both kind of fits 
provide us with compatible values of -Pg = °°(|<?| = 1>^)- 

We will use equation [12] to fit the values of P B =oc (q = 1,T) at different temper- 
atures. The argument we have done says nothing about the behavior of AP B (q = 

0, t,T), while APg(g, t,T) for q ^ follows the same scaling law valid for AP B (q = 

1, t,T). 

We have checked that our numerical data are well fitted by the behavior of equation 



12. We show two plots corresponding to two different temperatures: T = 0.7 and 
T = 0.35, in figure [^. In the plots we also draw our best fits to to the Ansatz [1^. 

The fit of figure [3] for T = 0.7, for example, looks reasonably reliable. The P B 
depends very slowly on T, and our numerical data have as support a range of times 
as large as the range we need to extrapolate through. On the contrary the fit of |3] 
for T = 0.35 (again, the one at our lowest T value, i.e. our most difficult fit) looks 
less happy even if one would not expect very strange things to happen: since at 
such low T is becoming very small our data vary on a small support, and the t = oo 
point looks far away. The very large error that appears in figure [| is an effect of this 
phenomenon. 

In figure |] we show the values of P B =oc (q = 1, T) extrapolated up to t = oo versus 
T. Because of how difficult it is to get precise results at low T values, the evidence 
we have established is not very strong: still, figure ^ clearly shows that our present 
low T results are completely compatible with the T = results. 

We end this section by performing a last consistency check. We compare the 
values of (q B ) with the values of the correlation function (in both cases extrapolated 
to infinite time). Using equation [| we can write equation [5] as: 

(q 2 B )(T) = 0.125 + 0.785 (7^(1, T) . (20) 

In table |l] we show the values of (q B )(T) obtained by direct measurements of q B and 
the value obtained by summing up correlation functions, according to the right hand 
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Figure 3: logPjg(l, t, T) versus t 22 < T ) for T = 0.7 (rhombs) and T = 0.35 (squares). 
The dashed line (dotted line) is for the best fit to the Ansatz [12] for T = 0.7 (T = 0.35 
respectively) . 



side of equation £CJ: the agreement of the two quantities is very reasonable at all 
T values. The values of (q%)(T) have been extrapolated at infinite time by using 
equation [14. 



4 Conclusions 

The results of our investigation are indeed positive. Physical observables computed 
at T = are compatible with values computed on large lattices thanks to an off- 
equilibrium dynamics, and extrapolated from measurements at low, finite T. Our 
observables of choice are interesting observables, since they do not only show that 
our system has a non-trivial P(q), but also that this non-trivial P(q) is the "bona 
fide" effect of replica symmetry breaking, and not the effect of some fancy kind of 
interfaces. 

We have been able to use very low T values, while still keeping a good control 
over t — > oo fits. When we lower T and eventually T — > the data are collapsing, 
in the reasonable accuracy given by the large statistical error, to the T = results. 
It is important to remark, also as far as planning future, more precise studies is 
concerned, that it is really very difficult to run accurate numerical simulations in the 
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Fi gure 4: Pg °°(1,T) versus T. The horizontal line corresponds to the confidence 
limit obtained in ground state computations. 



T 


(Q 2 B)(T) 


rhs of equation |20 


0.7 


0.585 ±0.002 


0.576 ±0.004 


0.6 


0.596 ±0.009 


0.594 ±0.004 


0.5 


0.57 ±0.02 


0.60 ±0.01 


0.35 


0.61 ±0.03 


0.64 ±0.02 



Table 1: Test of relation [H| 



region T < 0.3 with the nowadays computers (for example at T = 0.1 the dynamical 
critical exponents is of the order of 60). 
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